BMC Ecology and Evolution

Despite their currently scarce taxonomic diversity, rhinocerotoids were among the most evolutionarily successful clades of large herbivorous mammals throughout the Caenozoic. Within an array of morphological adaptations, the cranial horn represents the clades most remarkable and name-giving feature, imposing significant morpho-functional demands on the skull. In this paper, patterns and drivers of the rhinos cranial evolution were investigated for the first time by compiling a three-dimensional geometric morphometric dataset of living and extinct species. Multiple aspects of morphological evolution were explored using phylogenetic comparative methods, including craniofacial evolutionary allometry (CREA - tendency of larger species to have longer faces), phylogenetic constraints, and tempo of evolution. Cranial shape variation linked to horn presence and size evolved under a phylogenetically constrained framework, with a notable transition from hornless to horned species. Rhinos significantly deviated from CREA, likely due to diversity of cranial forms and proportions that evolved independently across the clade in response to varying horn morphologies and dietary habits. Hornless and horned rhinos exhibited similar rates of cranial evolution, and shape variation was obtained through multiple episodes of accelerated evolution. This highlights the role of morphological innovations and Caenozoic global cooling events in the emergence of phenotypic diversity.

The introduction of alien species into new habitats stands as a pressing economic and ecological challenge but it is also essential for unveiling evolutionary processes. The introduction of the Argentine ant (Linepithema humile) led to the spread of a single supercolony through different continents and thousands of kilometres like in Europe, from Northwest Spain to Greece. It was assumed that the high invasiveness of the species mainly relied on the lack of agonism among colonies, an effect derived from its introduction. However, recent studies suggest that local adaptations and evolutionary divergence could involve the disruption of the Argentine ant "empire" into a mosaic of competitive colonies. We investigated how isolation affects population divergence by comparing mainland and island populations in two distant regions colonized in Spain and in Greece with morphology, agonism, cuticular hydrocarbons, and genetic diversity of ant workers. Our results showed that all colonies sampled belonged to the most spread supercolony in Europe (main supercolony) except one sampled in Crete (Heraklion; Greece), which resulted to be a supercolony not registered in Europe. The Heraklian supercolony showed a different chemical and genetic profile and hostile agonism towards the other Greek colonies. Differences between islands and mainland colonies belonging to the main supercolony were higher in Galiza than in Greece. Surprisingly, the chemical profile of the Cretan colony belonging to the main supercolony showed more similarity with the Galizan colonies than with the Greek mainland, suggesting that L. humile may have been introduced into Greece through this island instead of the mainland. Our study suggests that local adaptations in Argentine ant colonies can trigger competition between colonies. Our data strongly support the existence of a candidate supercolony which highlights either ongoing introductions of L. humile in Europe or gaps in our understanding of its metapopulation dynamics.

BackgroundSynonymous mutations, once considered neutral, can affect translation efficiency through mRNA folding and splicing, generating codon usage bias. This bias is often linked to genomic GC content, which also influences gene regulation. In the parasitoid wasp Lysiphlebus fabarum, GC content was previously shown to shift between developmental stages, with larvae showing higher GC than adults. Whether this phenomenon is widespread among insects remains unknown. ResultsTranscriptomic data from six insect species spanning Diptera, Hymenoptera, and Lepidoptera was used to compare GC content between expressed genes in larvae and adults. In five species, larval transcripts exhibited higher GC content than adult transcripts. Differential expression analysis revealed that stage-biased genes displayed consistent GC shifts, and orthologous gene families with representatives across species showed particularly GC-rich larval-biased genes in Hymenoptera and Diptera. At the genome scale, modeling in 317 insect species demonstrated an association between parasitic lifestyle and reduced mean GC content in Hymenoptera and Diptera, providing a possible ecological explanation for AT-rich genomes. ConclusionsOur results show that GC content is dynamic across developmental stages, independent of overall genome composition. Stage-specific GC enrichment may reflect adaptive codon usage optimizing translation during energetically demanding life-history stages such as larval development. Furthermore, the association between parasitism and reduced genomic GC highlights how ecological lifestyle might with genome content and evolution. Lastly, this work identifies candidate genes underlying stage-specific GC bias and provides new insights into the interplay between molecular evolution, development, and parasitic adaptation in insects.

Hearing in Lepidoptera mediates a range of ecologically important behaviours, including mate communication, predator avoidance, and acoustic signalling. In moths, the evolution of predator-prey interactions with bats has further shaped hearing through a sensory arms race, with repeated co-option of auditory organs to detect and evade echolocating predators. Despite significant prior characterization of the neurophysiology and behaviour of hearing in moths, the genetic basis of hearing is poorly understood in most insects. In this study, we identify a core set of putative auditory genes in Lepidoptera using a combination of homology-based searches from Drosophila and evolutionary rate analyses. We find 56 genes present across all species and investigate whether gene copy number varies among non-hearing and hearing lineages and among 3 different ear types. We discovered seven genes associated with ear type and one with ear presence, but did not find significant losses in gene copy number in non-hearing species. We identified three genes (btv, Dnai2, and nompB) with strong evidence of selection in hearing clades and five genes with weaker evidence of selection. We discuss the potential roles of btv, nompB, and Dnai2 in ciliary transport and the aging of hair cells, as well as the possibility of actively amplified hearing. Our study serves as a primer and resource for further gene mining and functional testing of auditory genes in moths and other insects.

The between sex genetic correlation for traits has long been hypothesized as a constraint to the evolution of sexual dimorphism. Both empirical and theoretical work has suggested that this constraint is influenced by genotype-sex-environment interactions. We examine genotype-sex-environment interactions in both sexually exaggerated and non-exaggerated legs of Drosophila prolongata, to examine the role of organismal condition on evolvability of an extreme trait. We employed a nested full-sib half-sib crossing design, providing food either ad libitum, or restricting food during larval growth, to each brood. When provided food ad libitum, inter-sex genetic correlations between traits is high and positive, whereas under food restriction this correlation substantially weakens, with a modest negative sign. Similarly, comparisons of the G matrix across sexes becomes less associated under food restriction. We discuss these results in the context of the growing appreciation of the factors that facilitate sex-specific evolutionary change.

True crabs (Brachyura) are among the most iconic marine arthropods, representing noteworthy examples of morphological and ecological disparity. A striking feature of brachyurans are their anterior pincer-like appendages: chelipeds. These structures showcase a large diversity of morphologies that reflect ecology and overall multifunctionality. Yet, a comprehensive assessment of appendage functional morphology within phylogenetic and ecological trait contexts has never been attempted. By combining 3D geometric morphometrics, finite element analyses, multilocus molecular phylogeny, and ecological trait data for 80 crab species, including three fossil forms, we unveil a complex evolutionary history for crab chelipeds. Despite extreme shape diversity amongst chelipeds, stress distributions are very similar across taxa and hint a many-to-one pattern. High concentrations of chelipeds within constrained morphospace regions associated with peak pinch forces illustrates that brachyuran morphologies optimised for shell crushing may have arisen in the Cretaceous. Deviations from this morphospace highlight the diversification of non-shell-crushing life modes and the influence of sexual selection on appendages. Neither cheliped shape nor pinch force show phylogenetic signal. Together these results indicate that the evolution of cheliped shape is closely associated with, and inferred to have been strongly influenced by, crab ecology, biomechanical needs and sexual selection. SIGNIFICANCE STATEMENTChelipeds, the pincer-like claws of crabs, are among the most morphologically diverse appendages within Arthropoda, yet the evolutionary forces driving this diversity remain poorly understood. By integrating 3D geometric morphometrics, biomechanical modelling, molecular phylogeny, and ecological data across 80 crab species including fossil forms, we demonstrate that cheliped morphology is driven by ecology, biomechanical demands, and sexual selection rather than phylogenetic relatedness. The multifunctionality of these structures produces strong evidence for many-to-one mapping of form to function. Morphologies optimised for durophagy appear to have originated in the Cretaceous, with subsequent diversification into manipulative and sexually selected forms from a morphologically flexible foundation. These findings demonstrate that cheliped diversity reflects a complex interplay between ecological specialisation, biomechanical optimisation, and sexual selection across Brachyura.

Ecogeographic rules predict that endotherms inhabiting colder environments exhibit larger body sizes and shorter appendages to reduce heat loss. However, these hypotheses have largely been tested using size-based traits, leaving it unclear whether climatic adaptation also alters body proportions that more directly influence surface-area-to-volume ratio. Here, I examined ecogeographic variation in body shape, body size, and relative limb lengths across climatic gradients in three western North American mustelids: American ermine (Mustela richardsonii), American mink (Neogale vison), and Pacific marten (Martes caurina). I then modeled how variation in body shape and size affected surface-area-to-volume ratio. Body size increased in colder climates in all three species, consistently supporting Bergmanns rule. In contrast, only ermines became stouter in colder climates, whereas body shape remained unchanged in minks and martens. Patterns of appendage variation also differed: ermines and martens exhibited relatively longer limbs in colder climates, contrary to Allens rule, whereas only minks exhibited relatively shorter hindlimbs. Modeling showed that increased body size reduced surface-area-to- volume ratio in all species, whereas body shape contributed only in ermines. These findings demonstrate that body shape, body size, and appendage lengths respond independently to climatic gradients, with species-specific ecological and functional constraints shaping thermoregulatory adaptation.

Color polymorphism can facilitate local adaptation, maintain intraspecific diversity, or reflect early stages of speciation. Clownfishes (Amphiprion spp.) typically display a simple black, orange, and white pattern, but the saddleback clownfish (Amphiprion polymnus) shows striking variation in melanism and the number of vertical bars, which are thought to play a role in species recognition. In 2024, a revision on iNaturalist split A. polymnus into multiple species based solely on color pattern and geographic range. This raises the question of whether these morphs represent true species or intraspecific polymorphism, which we tested using genomic and image-based data. We sampled 97 individuals from seven populations across the species range and quantified color patterns from standardized photographs. Phenotypic and genomic analyses reveal a complex pattern of divergence. Image analysis identified three distinct phenotypic clusters, with A. polymnus, A. annamensis, and A. laticlavius each showing consistent differences in saddle shape and vertical bar extent. ADMIXTURE resolved three distinct genetic groups corresponding to the morphs. Pairwise FST (0.54-0.71) and dxy indicate extremely high differentiation between A. polymnus and A. annamensis, consistent with species-level divergence, whereas A. laticlavius shows much lower differentiation from A. polymnus (FST 0.09-0.18) and higher differentiation from A. annamensis (FST 0.64-0.66). Overall, phenotypic and genomic data show structured variation, but the status of A. laticlavius remains ambiguous. Our study reveals clear and structured divergence across the full range, yet the taxonomic interpretation of this variation remains inherently challenging. The key question remains: do these patterns reflect a single polymorphic species or a complex of closely related species?

BackgroundLyme borreliosis is a common tick-borne disease in Europe caused by spirochetes of the Borrelia burgdorferi sensu lato complex, including Borrelia afzelii, which is maintained in nature through interactions with rodent reservoir hosts. These spirochetes have evolved several surface proteins to manipulate rodent host immunity, some of which remain polymorphic in Borrelia populations. Among these proteins, OspE, which binds the host complement-regulating factor CFH to evade destruction by complement, is one of the most variable. Yet, what evolutionary forces maintain this polymorphism is not well understood. Motivated by a recent discovery of CFH polymorphism in the bank vole (Clethrionomys glareolus), the main reservoir host of B. afzelii, we hypothesized that the polymorphism is maintained by host-parasite coevolution involving specific associations between host and parasite genetic variants. MethodsWe analyzed associations between bank vole CFH alleles and B. afzelii OspE variants across three datasets sampled in Poland. Selection acting on OspE was evaluated using omegaMap. Host-pathogen genotype associations were tested using partial redundancy analysis (RDA), and co-structure was assessed using co-correspondence analysis (CoCA). ResultsWe found that OspE evolves under positive selection, however, we found no evidence for an association between OspE and host CFH variants at the individual level based on RDA or at the population level based on CoCA. ConclusionsDespite evidence of positive selection acting on OspE, we found no support for specific genetic matching between B. afzelii and its bank vole host at the CFH-OspE interface. These results suggest that the evolution of CFH and OspE may be shaped by broader selective pressures, potentially including interactions with multiple host species.

The origin of novel phenotypes can influence access to new ecological resources, which may have positive, neutral, or negative effects on subsequent phenotypic diversification. In this study, we tested the macroevolutionary consequences of a pair of putative functional innovations occurring in deep-sea fishes of the order Stomiiformes. Integrating phylogenetic comparative methods, micro-CT scans, and external body measurements, we recover a mosaic of diversification trends associated with these innovations. We found some evidence for elevated evolutionary rates in tooth morphology associated with the predatory dragonfishes, which possess a gap between their vertebral column and skull that exposes the notochord and enables neck-like flexibility. However, a second novelty building upon the first, a functional neck joint enabling extreme cranial kinesis, was linked to faster rates of skull evolution. Our results suggest that innovations that help shift ecological roles and overcome functional constraints related to those roles, like gape-limitation in prey depauperate habitats, may play an important role in promoting phenotypic diversification. This work builds on a growing body of evidence highlighting how the deep sea promotes phenotypic diversity, generating the extreme forms that are celebrated by scientists and the public alike.

Heatwaves can affect species abundances by changing how species interact with one another in local communities. These effects can be complex and remain poorly understood, especially in cases where the outcome of species interactions is contingent on the history of species arrival. We studied how heatwaves affect interactions between the yeast Metschnikowia reukaufii and the bacterium Acinetobacter nectaris, both commonly found in the floral nectar of Diplacus aurantiacus, a hummingbird-pollinated shrub native to California. The microbes were introduced to artificial nectar in different orders of arrival in the presence or absence of simulated two-day heatwaves. We found that heatwaves made yeast-bacterium interactions more contingent on arrival history, thereby causing large variation in nectar acidity, a factor known to affect hummingbird preference and seed production. In the absence of heatwaves, Acinetobacter always became abundant regardless of arrival history, suppressing Metschnikowia and reducing nectar pH. In contrast, in the presence of heatwaves, Acinetobacter dominance depended on arrival history and heatwave timing. If Acinetobacter arrived after Metschnikowia during a heatwave, Metschnikowia suppressed Acinetobacter substantially enough to keep nectar pH at a high level. These results suggest that heatwaves shift species interactions from determinism to historical contingency, with both taxonomic and functional consequences.

The study of historical documents from the original excavation led by J. M. Cruxent, of museum collections, as well as results of fieldwork, are combined to provide a revised assessment of the archeologically site of Cucuruchu, which together with those of Taima-Taima and Muaco, document megafaunal lithic remains associated from the Late Pleistocene in northern Venezuela. Six vertebrate taxa are reported for Cucuruchu: five megaherbivores (a giant sloth, Eremotherium laurillardi; a glyptodontid, Glyptotherium cf. G. cylindricum; a macrauchenid, cf. Xenorhinotherium bahiense; a proboscidean, Notiomastodon platensis; and an equid, Equus sp.), and a turtle (Kinosternon sp.). Two fragments of El Jobo projectiles were found in situ lying among megaherbivores bones, while another unidentified projectile fragment was also found in the same layer. Although no direct evidence of human exploitation on these megaherbivores remains has been found, the presence of El Jobo projectiles in the undisturbed fossil-bearing layer some of which are associated with fossil remains, could be a suggestive of such interaction. Two new radiocarbon dates indicate ages of at least ca. 15.3-16.0 kyr cal BP for the fossil-bearing layer, contrasting with the mid to late Holocene age originally reported. Cucuruchu is a window to the past and supports the presence of hunter-gatherers of the El Jobo culture in the region at the end of the Pleistocene.

The nautilus (Nautilus sp.) is an early-branching cephalopods. It retains several conchiferan synapomorphies, including an external planispiral biomineralized shell. The shells complex structure allows it to withstand hydrostatic pressure, control buoyancy, and protect against external hazards. In this study, we comprehensively examined shell microstructures across different shell components and regions representing various ontogenetic stages in two adult museum shell specimens. We found that the nautilus shell is composed of five microstructural types (spherulitic, prismatic, nacreous, semi-prismatic, and irregularly oriented prismatic structures) organized into layered architectures within individual shell components and coordinated across the shell as an integrated system. Our observations highlight transitions between distinct microstructures within and across shell components and local variation within individual components such as the dorsal and ventral shell walls, suggesting that these patterns may contribute to shell strength and overall mechanical performance. Variation in caecum morphology suggests that this structure may be developmentally plastic and subject to relatively relaxed structural constraints. These findings show that the Nautilus shell is an integrated biomineral system in which diverse microstructures are organized across shell components to meet functional demands and provide the mechanical strength needed for survival.

High investment in current reproduction can limit future reproductive opportunities, but how selection shapes these hormone-mediated traits remains poorly understood. Androgens can mediate male reproductive investment, and in three-spined stickleback (Gasterosteus aculeatus), exert antagonistic effects on breeding effort versus spermatogenesis. To understand how shifts in reproductive strategy shape this tradeoff, we compared testes transcriptomes and androgen production between two recently diverged stickleback ecotypes that differ in reproductive strategy: the ancestral "common" ecotype, which provides paternal care, and the non-parental "white" ecotype, which has lost paternal care and prioritizes mating effort. During typical breeding, testes gene expression differed little between ecotypes. However, under prolonged summer-like conditions, testes gene expression diverged substantially. Common-biased genes were enriched for meiotic functions and spermatogenic cell type markers, suggesting commons had initiated spermatogenesis while whites had not. Instead, whites expressed higher levels of steroidogenic candidate genes and released significantly more 11-ketotestosterone than commons, indicating sustained investment in current reproduction. F1 hybrids released 11-ketosterone at intermediate rates, suggesting a genetic basis for this divergence. Sustained androgen production in whites may possibly delay the transition into spermatogenesis, limiting investment in future reproduction. These results illustrate how selection on hormonally-integrated traits can drive rapid divergence in life history strategy.

Late Pleistocene sea-level fluctuations and glacial corridors intermittently connected the island of Ireland to Britain and continental Eurasia, shaping patterns of megafaunal dispersal and occupation. However, the scarcity of genetic data from ancient Irish fauna has limited our understanding of their demographic histories and relationships to continental populations. To address this, we generated sedimentary ancient DNA sequences from Castlepook Cave in southwest Ireland, detecting twelve ancient taxa, including two without zooarchaeological records at the site. We recovered the first mitochondrial sequences from Irish cave hyenas and woolly mammoths, providing new insights into their maternal population history: cave hyenas carried mtDNA haplogroup A1, previously identified in Late Pleistocene European populations, while woolly mammoths belonged to clade III/B2, which was replaced in Europe at approximately the same time. We also identify two mitochondrial clades (1b and 2), consistent with pre-Last Glacial Maximum mitochondrial lineage turnover in Ireland. Together, these results indicate that Irish megafauna formed a biogeographical continuum with their conspecifics in Europe, extend the known geographical range of several mtDNA lineages, and support faunal connectivity at the northwestern edge of Europe prior to the Last Glacial Maximum.

Mutation is the ultimate mechanism that produces genetic novelty, and thus a central ingredient of evolution. Mutation rates are therefore thought to be tuned by natural selection, for example to optimize a delicate balance between the generation of adaptive diversity and the accumulation of deleterious mutations. As this selection occurs over very long time scales, models and simulations have been powerful tools to understand how mutation rate evolves and which factors influence it. Most simulation methods are nevertheless limited by the over-simplicity of the genotype-to-phenotype map they feature, especially regarding the encoding of mutation rate. We modified Aevol, an evolutionary simulator inspired by bacterial genomics with a realistic genome structure and a complex genotype-to-phenotype layer, to allow organisms to evolve genes coding for higher replication fidelity. This setup permits several degrees of realism absent in other models: mutation-rate modifier genes themselves experience a realistic distribution of effects of mutations and diminishing- returns epistasis, similarly to fitness modifiers. Moreover, a lower mutation rate comes with the trade-off of a larger genome to encode the genes improving replication fidelity. We use this setup to test hypotheses regarding the evolution of prokaryotic mutation rate, and its link with genome size and genetic drift. We found that evolution systematically increases replication fidelity, even when this results in lower fitness. We highlight two factors which limit the mutation rate decrease: genetic drift and the supply of gain-of-fidelity mutations.

This research addresses a central question in Palaeolithic research: how hunter-gatherer mobility was structured across space and time in the Cantabrian Region (northern Iberia), which has human occupation evidence spanning from the Middle Pleistocene through the Holocene. A multidisciplinary framework integrating primarily {delta}3S isotope values, combined with {delta}{superscript 1}3C and {delta}{superscript 1}N, palaeoproteomics, Bayesian age modelling, palaeoclimatic reconstruction, isoscape mapping, and ecological diversity was developed. A total of 905 animal bone collagen samples, with evidence of anthropogenic modifications, from 16 key archaeological sites from the Mousterian to Mesolithic (Marine Isotopic Stage 5 to 1, between 100-7 ka BP) were analysed, permitting the reconstruction of spatial patterns of resource exploitation and human mobility. The {delta}3S isotope values show weak, inconsistent relationships with climatic proxies, suggesting that sulfur signatures are primarily driven by geographic and ecological factors rather than climate. Strong spatial trends are observed, with higher {delta}3S values in coastal zones and lower values inland. Diachronic trends reveal marked shifts in human mobility: smaller ranges during the Mousterian, increasing mobility through the Chatelperronian and especially the Aurignacian, followed by reduced mobility in the Gravettian and Solutrean, and renewed territorial expansion during the Magdalenian and, likely, the Azilian. In contrast, the Mesolithic is characterised by decreased mobility and thus increased territoriality in both coastal and inland contexts. Faunal isotope values and isoscape predictions reveal that some animals were acquired beyond local foraging ranges during the Palaeolithic, particularly in inland regions with lower {delta}3S values. Isotopic niche analyses indicate partial interspecific overlap consistent with ecological flexibility. Macromammal and micromammal diversity exhibit contrasting patterns, with a significant negative correlation in Simpson and Shannon indices. Macromammal diversity correlates negatively with {delta}3S values, linking increased hunting diversity to expanded catchment areas and longer-distance foraging, whereas micromammal diversity shows positive correlations with {delta}3S, {delta}{superscript 1}3C and {delta}{superscript 1}N reflecting stronger climatic influence. Overall, these results demonstrate that hunter-gatherer behaviour in northern Iberia during the Middle and Late Palaeolithic was highly dynamic, combining logistical and residential strategies that shifted in response to changing environmental conditions, resource distributions and cultural adaptations.

The unprecedented rate of climate warming threatens many species, and assessing their vulnerability to climate change represents a critical challenge in conservation biology. The Apollo butterfly, an emblematic mountain species, is expected to be impacted by climate change. Here, we analysed thousands of SNPs from 101 localities across Apollo French distribution. We identified 93 SNPs strongly associated with climate variation using five genotype-environment association analyses. We forecasted future climate maladaptation of French Apollo populations using four genomic offset methods and integrated these results with neutral and adaptive genetic diversity, genetic structure and adaptive climatic niches to infer their vulnerability to climate change. Jura and Alps populations exhibited the lowest risk of vulnerability to climate change, with low genomic offsets, high genetic diversity and connectivity, whereas Auvergne populations showed the highest genomic offsets and lowest neutral and adaptive genetic diversity. Only a reduced percentage (<1%) of the current distribution is predicted to face climatic conditions outside the current range, suggesting that adaptive variability required to adapt to future climates may already be present, and that assisted gene flow could represent an effective conservation strategy. Finally, we discuss some of the main challenges of genomic forecasts, particularly for declining non-model species.

Sharks exhibit extraordinary morphological diversity across a wide range of ecological niches, yet large-scale, high-resolution digital datasets of their internal anatomy remain limited. Here we present an open-access 3D shark anatomical repository derived from published X-ray computed tomography (CT) data, featuring manually segmented and systematically annotated models of the chondrocranium, visceral arches, axial skeleton, musculature, and viscera in standard STL format. The dataset comprises 117 individuals, representing 72 species across 25 families and all nine extant shark orders, with 115 full-body reconstructions and two head-only models. This open-access dataset offers a comprehensive resource for comparative anatomy, biomechanical simulations, evolutionary developmental biology and biomimetics research of extant sharks.

Yeasts thriving in floral nectar have a proven ecological engineering capacity, as they modify the internal environment and the rewards of individual flowers in ways that influence how plants and pollinators interact. This floricentric perpective, however, produces only a partial view of the potential role of yeasts in plant-pollinator interactions, and little is known on their possible engineering effects on whole inflorescences. This study assessed yeast ability to shape features of the vertical inflorescences of Gladiolus illyricus (Iridaceae), a Darwins syndrome plant (acropetalous protandrous flowers pollinated by insects foraging from bottom to top). When yeasts were excluded from inflorescences a declining gradient in nectar sugar concentration built up from the oldest, female-stage flowers at bottom up to the youngest, male-stage ones at top. This plant-intrinsic gradient was reversed in inflorescences with yeasts, where sugar concentration increased from the female-stage flowers harboring dense yeasts population to male-stage ones with few yeasts. The two species of yeasts involved (Metschnikowi reukaufii and M. gruessii) differed in their effects. Results add support to the growing consensus that specialized nectar yeasts can obfuscate intrinsic patterns of intra and interspecific variation in flowering features of plants at the level of individual flower, inflorescence, population and plant community.